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ABSTRACT 

Using partial wave unitarity and the observed density of the Universe, we 
show that a stable elementary particle which was once in thermal equilibrium 
cannot have a mass greater than 340 TeV. An extended object which was once 
in thermal equilibrium cannot have a radius less than 7.5 x 10~ 7 fm. A lower 
limit to the relic abundance of such particles is also found. 
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The idea that the dark matter known to exist in galactic halos consists of 
some, as yet undiscovered, stable massive particle has received a great deal of 
attention in recent years. Dozens of particle candidates have been suggested and 
new ones are constantly being proposed. Most of these dark matter candidates 
have relic abundances which are calculated in the “Lee- Weinberg” manner and 
model parameters are typically adjusted to allow their density today to be near 
critical density, fix ~ 1. In this Letter we wish to point out, that for almost any 
such particle which was once in thermal equilibrium and has an abundance deter- 
mined in this way, partial wave unitarity of the S matrix bounds the annihilation 
cross section in the early universe, which in turn bounds the relic abundance and 
the mass of the particle. In general we find that stable elementary particles with 
masses greater than around 340 TeV are very likely excluded. Extended objects 
with radii less than 7.5 x 10 -7 fm are also very likely excluded. 
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As an application of these limits, we note that the claim of Enqvist, tt al ., 
that there is no upper limit from cosmology on the mass of a stable Dirac neutrino 
cannot be true. While the mass upper limit we find is not rigorous and rather 
high, we still feel it may be of some interest because of its general nature. 

The relic abundance of a particle species, X, which was once in thermal equi- 
librium, is determined by its total thermally averaged annihilation cross section 
(<r(XA — > aU)v re i) at freeze-out. At high temperatures the number density of A’’s 
is roughly the same as the number density of photons, but as the temperature 
drops below the mass of the A, their number density drops exponentially. This 
continues until the total annihilation cross section is no longer large enough to 

maintain equilibrium and the A number density then freezes-out- - The number 
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density today is given roughly by 




1.07 x 10 9 (n + l)x" +1 GeV 1 3 x 10 2 ‘cm 3 /sec 


gV 2 m p i{av ie \)' f 


(<™rel) / 


( 1 ) 


where Dy = px/pcrit is the present average density of A’s divided by the critical 
density, \ < h < 1 is the Hubble constant in units of 100 km/sec/Mpc, Xf = 
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m x/Tf, Tf is the freeze-out temperature, g, % 107 is the effective number of 
degrees of freedom at Tf , and m p i = 1.22 x 10 19 GeV. 

Since the XX annihilations at freeze-out occur at non-relativistic velocities 
(u | -C 1), one can expand the cross section in powers of v 2 = v 2 e ^/4 and keep 
only the first (or first two) terms. In thermal averaging one replaces (u^j) = 6/xy 
and so in eq. (1) the cross section is written (crt; re i) = (< 7 i> rel )'x _n , where n pa- 
rameterizes the dependence of the cross section on x. The freeze-out temperature 

* , 3 

is given roughly by 

xy = In I? — (n + — ) In In B, (2) 

where B = .038gm p im x {av te \)' / y/g* and g is the number of degrees of freedom of 
the A particle. Typically xy ss 25 corresponding to u 2 e] /4 « 1/16 at freeze-out. 
Please note that t> re ] is not really a velocity, but is related to the flux factor. It 
is defined as u re i = 2u, and so 0 < ^el/4 < I- 

From eq. (1) we see that if (crv re j)y «3x 10~ 27 cm 3 /sec, then Q x h 2 > 1, 
which would be inconsistent with the “observation”, f2tot^ < 1, obtained from 
the age of the universe. Any particle model which predicts an annihilation cross 
section smaller than this critical value at v “ el /4 « ^ is therefore inconsistent with 
cosmology. We will now’ show that partial wave unitarity provides a maximum 
possible cross section and therefore a minimum possible Q x h 2 . Extremely mas- 
sive elementary particles and very small extended objects violate these bounds 
and therefore are inconsistent with cosmology. 

Consider the process a + b — * c -f d and the scattering matrix 

</|S|>) = m + i(2ir) 4 «V/ - P.) {/mo (3) 

where P t = p a + Pb and Pf = p c + pj. The T matrix can be expanded in partial 
waves using the helicity formalism 5 

(A c Arf|T(s, n)\X a X b ) = 87rs 1 /2 e ^(A-V) £( 2 J + l)dl x .(6) {X c X d \Tj(s)\X a X b ) , (4) 

J 

where A a , • ■ • , Aj are the helicities of particle a,- ■■ ,d. X = X a - X b , X' — X c — X^, 
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s is the Mandelstam variable, ft = (0, <j>) is the center-of-mass scattering angle, 
and d J xx , are the Wigner functions. 

Using matrix notation 6 (A c A,ijTj(s)|A 0 Aj ) ) = ( Tj)if and pk = diag(pi,p2i • ‘ ')•> 

where pk is the center of mass three-momentum of particle system i, /, etc. partial 

, , 6 

wave unitarity of the S matrix can be written 


Tj - Tj = 2iTjpTj. (5) 

Defining Sj = 1 4- 2ip l / 2 Tjp 1 / 2 , we see that partial wave unitarity can also be 
written SjSj = 1 or 

\Sel,j \ 2 + £ ( 6 ) 

f 

where S ( i j stands for the elastic channel, i = /. The next step is to define 
S € \ j = where 8j is a real phase shift and tjj is an inelasticity factor, 

0 < T)J < 1. Then |5 e/ ,j| 2 = tj], and E/l 5 *¥/,j| 2 ■* 1 ~ Fina11 ^ usin S 
T e i,j = (Seij-l)/(2ip) and T/^j = Sf^i,j/{2i y /pipJ), and the standard formula 
for the unpolarized cross section in terms of partial waves a = where 


vj = 


4tt(2 J + 1) 
(2s a + l)(2sj, -f 1) 



(7) 


we 


find the result of Pilkuhn 
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0r,J = 4 tt 


(2J + 1) 

(2s a + 1)(2 s 6 + 1) 



*(2 J + i)(i - n]) 
Pi 


( 3 ) 


Here cr r j is the u reaction” cross section, that is, the total cross section minus the 
elastic piece. It has a maximum when rjj = 0, so we conclude that 


a j(a + b — *■ c + d) < 


7r(2J + l) 


P? 


(9) 
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In the eaxly Universe, 


Pi 


= E l 


2 

m y 


m X v i\ 




4 (! - v 2 x J4) 


so crjv K \ < ( <7j) max v re i , where 


(&J )mai^rel ^ 


4?r(2J + 1) 

m\-v re i 


3 x 10 22 (2J + l)cm 3 /sec 
(m A -/TeV) 2 


( 10 ) 


In order to apply the limits of eq. (10) to the annihilation in the early Uni- 
verse we need to determine which partial waves contribute. After summing over 
helicities, the angular dependence, cos 6, which indicates the partial wave, enters 
the cross section only through the Mandelstam variable 


t -m\ + m 2 - 2 E a E c + 2 p c p a cos 6 
=m 2 + m 2 - 2 E a E c 4- 2 p c cos 0 m A t > re \/2 + 0(u^ el /4). 


( 11 ) 


So there is a factor of u re j appearing with every factor of cos 6. In the expansion 
of the annihilation cross section in powers of 4,/4 1/16, the lowest order 

term 0((u 2 e |/4)°), therefore has no angular dependence and must be a J = 0 
partial wave. The J = 1 partial wave is smaller by a factor e 2 e) /4, and the 
higher partial waves are further suppressed. In fact, since partial wave unitarity 
must hold for any value of i> 2 e) /4, and when u 2 el /4 increases, the maximum cross 
section, eq. (10), decreases, the J = 0 bound, taken when 4l/4 1/16, is not 

as stringent as possible. The J = 1 maximum cross section also decreases for 
larger t?^ el /4, and more importantly, the term in the actual cross section of order 
v rel/4> increases. If the J = 1 bound is satisfied for a larger value of t£l/4, for 
instance <4/4 1 /2, then the J = 1 partial wave is below the bound by a factor 

of 8 -3 / 2 .04 by freeze-out. We conclude that it is more than adequate to use 

only the J — 0 partial wave in finding a bound. 
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Now we use eqs. (1), (2), and (10) to bound Q x h 2 and m x . Including only 
the n = 0 part of the cross section and replacing u re j = -y/6/xy, we find that 

n x h 2 > 1.7 x 10~ 6 ^x7(m x /TeV) 2 (12) 

for a Majorana fermion with g = 2. For a Dirac fermion, Q x h 2 is a factor of two 
larger. Now using Q x h 2 > 1, we find the mass limit 

m x < 340 TeV, (13) 


and x f « 28. Eq. (13) was found for a Majorana fermion. The limit for a scalar 
particle is similar, while for a Dirac fermion is about a factor of >/2 smaller, that 
is, m x < 240 TeV. This is the main result of this Letter. 


Another, more conservative, way of finding the mass bound is to assume 
that the cross section, eq. (10), holds throughout the period of annihilation and 
freeze-out. In this case, the t' re ' factor affects the thermal averaging and the 
integration from freeze-out to today. The thermally averaged maximum cross 
section becomes 


(( & J ) max v rel) 


4tt(2J + 1) 


m 


X 



(14) 


and the relic abundance is given by eq. (1), with n = -1/2, 


Sl x h 2 > 


6.0 x 10~ 7 

(2 J + 1) 


z} /2 (m*/TeV) 2 . 


(15) 


The freeze-out temperature is the same as before with (crv Te \)' multiplied by a 
factor of (6/7I-) 1 / 2 . (We set n = -1/2 in eq. (2), both now and before, since 
the xj/” here is just an algebraic factor.) Using these formulas, the mass limit 
becomes m x < 550 TeV. This is probably an overly conservative bound since 
one does not expect av K \ oc tQ 1 for annihilation channels in a nonrelativistic 
expansion. 
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However, we do not claim that the derivation leading to eq. (13) is rigorous, 
or that exceptions cannot occur. For example, elastic scattering via f-channel 
exchange of a massless particle gives rise to a term in the matrix element pro- 
portional to t~ l oc (1 — cos#) -1 . Naively expanding this would suggest that 
all partial waves contribute to the term of lowest order in u 2 el /4. The problem, 
in this case, is that we are outside the Lehmann ellipse of convergence, and the 
partial wave expansion not valid. Fortunately, in annihilation, the mass of the 
annihilation product must be less than mx , and the partial wave expansion con- 
verges, giving nicely the results we claim above. Another possible exception, 
which we do not consider very likely, is that the coefficients of the partial wave 
expansion contain factors of ( s — 4m^) -1 oc v~J, in just such a way as to cancel 
the v 2 e , factors associated with the cos 2 9 factors. For elastic scattering, it can 
be proved that this cannot occur (Ref. 6, page 291), but we have been unable to 
complete the proof for the inelastic case. This may be related to the possibility 
of s-channel poles, which can cause another possible exception to our limit. A 
factor of (s — m 2 ) -1 , with m, = 2m x will give an additional factor of t> -2 , in 
which case partial waves up to J — 2 need to be included in our maximum cross 
section, and the mass limit weakens. However, we feel that such a pole is unlikely. 
It requires not only an exchange particle of precisely twice the mass of the .V, 
but also that the exchange particle be nearly stable. The width of the exchange 
particle will dominate the pole unless it is very small, and since the exchanged 
particle is more massive than the X , and has decay channels into lighter particles, 
we consider this possibility remote. 

We note that the mass limit, eq. (13), involves a mass somewhat higher than 
typically considered in particle dark matter model building. But since the bound 
is so general we feel it may be of some use. As an example, we can immediately 
apply it to candidates which appear in the literature, such as the Dirac neutrino. 

The Dirac neutrino was the first dark matter particle considered and very 
early Zeldovich claimed a range of neutrino masses, 3 GeV < m v < 3 TeV, as 
being cosmologically acceptable. His upper bound was based on neutrino annihi- 
lation into fermions through Z boson exchange. This cross section is proportional 

. 2 

to m -2 in the high mass limit. However, Enqvist, Kainulainen and Maalampi 
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noted that the W + W~ channel, among others, open up for very massive neutri- 
nos, and that these new channels dominate the cross section in the high mass 
limit. In fact, they claimed that because the cross section keeps growing as m v 
increases, there is no upper limit from cosmology on Dirac neutrino masses. This 
claim is clearly contradicted by the bound for Dirac fermions given just after 
eq. (13). Yet we do not believe that that the cross section of Enqvist, et al. is 
in error. We believe that the solution to this puzzle is that in the Standard 
Model, where neutrinos get their mass by the Higgs mechanism, as m„ — ► oo, the 
neutrino Yukawa coupling becomes large and perturbation theory breaks down. 
Another way of saying this is that the higher loop corrections become important 
in this limit and the tree level calculation of Enqvist, et al. is not applicable. In 
fact, by using unitarity to bound the largest eigenvalue of the scattering matrix, 
Chanowitz, Furman and Hinchliffe, showed that the breakdown of perturbation 
theory occurs at around m„ «s 1 TeV, far below the limit we set. The breakdown 
of perturbation theory suggests that the neutrino becomes “strongly interacting” 
and could not exist as a free, stable state. In this case, the annihilation cross sec- 
tion would be governed by different physics, if the theory made sense at all. If, on 
the other hand, the neutrino for some reason stays “elementary”, we argue that 
our limit applies, giving an upper limit on the neutrino mass from cosmology, 
just as Zeldovich originally suggested (though at a different value). 

Finally, we should comment on the applicability of these bounds to extended 
objects. For these objects, higher partial waves will generally contribute to the 
nonrelativistic cross section, and the cosmological mass bound, eq. (13) does 
not apply; however, partial wave unitarity may still be used to limit the total 
annihilation cross section, and cosmology provides a constraint on the size of such 
objects. Consider an extended object with spin 0 and radius Rx • The highest 
partial wave that can contribute to the particle-antiparticle collision is roughly 
Jmax = 2mxv r eif2xi resulting in a maximum total cross section, 

a J max 

W W — £< 2J + 1) w lerfvtw, (16) 

m X V re\ J=rQ 

four times the geometric cross section. Using eqs. (1), (2), and (16), we can now 
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bound Qxh 2 and Rx ■ We find that 


> 


4 x 10 -15 x| 
(Ajr/fi n ) 2 


which leads to the bound 


■Rx > 7.5 x 1CT 7 fm. 


(17) 


(18) 


Here we used xj = 27 which was obtained from eq. (2) using mx = 1000 TeV; 
the radius limit, eq. (18), varies only logarithmically with mx- The limit for spin 
^ particles is more stringent by a factor of %/2. 

We point out that eq. (16) is valid only if J max ^ 1- • On the other hand, 
if Jmax "C 1, the cross section is bound by eq. (9) with J=0. Since freeze-out 
occurs when u te i ~ eq. (18) is reliable only when Rx l/mx, while an object 
with Rx <C 1 /mx must be considered point-like and its mass limited by eq. (13). 
Furthermore, we note that there is no major discontinuity in the overlap region, 
Rx rsj 1/mx, since the mass limit for point-like particles, eq. (13), is very nearly 
that which we would have obtained from the radius limit, eq. (18), had we used 
the Compton wavelength of the particle for Rx- 

Of course, if some process such as a quark-hadron or electroweak phase tran- 
sition, out-of-equilibrium decay of a massive particle, or inflation produces a 
significant amount of entropy after freeze-out, the relic abundance is diluted and 
our limits are weakened accordingly. Nevertheless, although our derivation is not 
rigorous, and exceptions may exist, we believe that the limit on mass, eq. (13), 
radius, eq. (18), and relic abundance, eq. (12) is of great interest and applies to 
many (if not most) dark matter candidates. 
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